Parallels between the dynamics at the noise-perturbed onset of chaos in logistic maps 

and the dynamics of glass formation 
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We develop the characterization of the dynamics at the noise-perturbed edge of chaos in logistic 
maps in terms of the quantities normally used to describe glassy properties in structural glass 
formers. Following the recognition [Phys. Lett. A 328, 467 (2004)] that the dynamics at this 
critical attractor exhibits analogies with that observed in thermal systems close to vitrification, 
we determine the modifications that take place with decreasing noise amplitude in ensemble and 
time averaged correlations and in diffusivity. We corroborate explicitly the occurrence of two-step 
relaxation, aging with its characteristic scaling property, and subdiffusion and arrest for this system. 
We also discuss features that appear to be specific of the map. 
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I. INTRODUCTION 

The erratic motion of a Brownian particle is usually 
described by the Langevin theory pj . As it is well known, 
this method finds a way round the detailed consideration 
of many degrees of freedom by representing via a noise 
source the effect of collisions with molecules in the fluid 
in which the particle moves. The approach to thermal 
equilibrium is produced by random forces, and these are 
sufficient to determine dynamical correlations, diffusion, 
and a basic form for the fluctuation-dissipation theorem 



In the same spirit, attractors of nonlinear low- 
dimensional maps under the effect of external noise can 
be used to model states in systems with many degrees of 
freedom. In a one-dimensional map with only one con- 
trol parameter /i the consideration of external noise could 
be thought to represent the effect of many other systems 
coupled to it, like in the so-called coupled map lattices 
0. Notice that the general map formula 



H+l 



x t + h M (x t ) + oft, 



(1) 



is a discrete form for a Langevin equation with nonlinear 
'friction force' term h^. In Eq. t = 0, 1, . . . is the 
iteration time, £t is a Gaussian white noise ((£t£e) = 
S t y), and a measures the noise intensity. 

An interesting option offered by Eq. is the study 
of singular states known to exhibit anomalous dynam- 
ics. For instance, the so-called onset of chaos in logistic 
maps is a critical attractor with nonergodic and non- 
mixing phase-space properties, and the perturbation of 
this attractor with noise transforms its trajectories into 
genuinely chaotic ones with regular ergodic and mixing 
properties. A case in point here would be to obtain the 



atypical dynamics near an ergodic to nonergodic transi- 
tion and compare it, for example, to that in supercooled 
liquids close to glass formation, where also an ergodic to 
nonergodic transition is believed to occur. The dynam- 
ical properties for such glassy states appear associated 
to two-time correlations with loss of time translation in- 
variance (TTI) and aging scaling properties, as well as 
to subdiffusion and arrest @, Q . The specific question 
we would like to address is whether there are properties 
shared, and if so, to what extent, by critical attractors in 
nonlinear low-dimensional maps and nonergodic states in 
systems with many degrees of freedom. 

Recently |||, it has been realized that the dynamics 
at the noise-perturbed edge of chaos in logistic maps 
shows similarities with that observed in supercooled liq- 
uids close to vitrification. Three major features of glassy 
dynamics in structural glass formers, two-step relaxation, 
aging, and a relationship between relaxation time and 
configurational entropy, were shown to be displayed by 
the properties of orbits with vanishing Lyapunov coeffi- 
cient. Interestingly, the previously known properties in 
control-parameter space of the noise-induced bifurcation 
gap (See Fig. ^| 6, ,7| play a central role in determining 
the characteristics of dynamical relaxation at the chaos 
threshold, and this was exploited to uncover the analogy 
between the dynamical and the thermal systems in which 
the noise amplitude a plays a role equivalent to a tem- 
perature difference T — To, where T g is the so-called glass 
transition temperature 0, • 

In Ref. p| only the properties of single-time functions 
(i.e. trajectories) were discussed and here we focus the 
analysis on two-time correlations. This would allow for a 
closer examination of the analogy with glassy dynamics 
as the study of the latter often centers on two-time cor- 
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relations. However, we also look at diffusion properties 
via the determination of the mean square displacement 
of trajectories in a suitably space-extended map. To lay 
emphasis on the similarities in the dynamics between the 
two types of systems we analyze the attainment of TTI in 
the correlations caused by the action of noise. We recall 
0, |j| that exposure of attractors to noise have the fea- 
tures of an activated process a mechanism that is usually 
considered in the interpretation of relaxation processes 
in glass formers. Also we make up a simple 'landscape' 
picture for the properties of the noise-perturbed map in 
order to compare it with that obtained from the mul- 
tidimensional energy landscape of supercooled liquids. 
Aware of our crude attempt to contrast the dynamics 
of a single map (although equivalent to a coupled ar- 
ray of such maps when described via a Langevin-type 
equation) with that of thermal systems we point out the 
main differences encountered. Specifically, the dynam- 
ics at the onset of chaos displays regular patterns absent 
in the known (experimental or computed) dynamics of 
molecular systems. These differences reflect the pecu- 
liarities of the period-doubling route to chaos displayed 
by unimodal maps. 

The structure of the body of the article is as follows: 
In Section 2 we recall essential properties of the logistic 
map under additive noise, e.g. the bifurcation gap and its 
time-dependent manifestation at the chaos threshold. In 
Section 3 we present results on ensemble-averaged two- 
time correlations for a > 0, the attainment of TTI, and 
the development of two-step relaxation as a — > 0. In Sec- 
tion 4 we focus on time-averaged two-time correlations 
and the occurrence of the characteristic aging scaling at 
the onset of chaos with a = 0. In Section 5 we present a 
repeated-cell map for diffusion at the onset of chaos and 
show the evolution in behavior from diffusive to subdif- 
fusive and finally to localization as a — ► 0. In Section 6 
we make final remarks. 



II. NOISE-PERTURBED ONSET OF CHAOS 
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We describe now the effect of additive noise in the dy- 
namics at the onset of chaos in the logistic map 

xt+i = 1 - \ix\ + aiu "I < x t < 1, < n < 2, (2) 

where in Eq. Q has taken the form (x) = 1 — x — 
fix 2 . For a = at the onset of chaos at /i c = 1.40115... 
the orbit with attractor initial condition xq = consists 
of positions arranged as nested power-law time subse- 
quences that asymptotically reproduce the full period- 
doubling cascade that occurs for (j, < fx c 0, hll. See 
empty circles in Fig. [21 This orbit is the last (the ac- 
cumulation point) of the so-called 'superstable' orbits of 
period 2™ which occur at \i = ~p n < fj, c , n = 1, 2, a su- 
perstable orbit of period 2°°. Superstable orbits include 
x = as one of their positions and their Lyapunov ex- 
ponent Ai diverges to — oo [||. A transient dynamics at 
the onset of chaos is observed for trajectories with initial 



FIG. 1: (a) Logistic map attractor. (b) Magnification of the 
box in (a), (c) Noise-induced bifurcation gap in the magnified 
box. 



position outside the Feigenbaum attractor but we shall 
not consider this in what follows. For a > the noise 
fluctuations wipe the fine structure of the periodic attrac- 
tors as the iterate visits positions within a set of bands 
or segments similar to those in the chaotic attractors, 
nevertheless there remains a well-defined transition to 
chaos at U c (c) where the Lyapunov exponent Ai changes 
sign 0, . The period doubling of bands ends at a fi- 
nite maximum period 2 N ^ as \i — * Mc(c) (see Fig. EJ:) 
and then decreases at the other side of the transition. 
This effect displays scaling features and is referred to 
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FIG. 2: Absolute values of positions in logarithmic scales of 
iterations t for various trajectories at the onset of chaos A J c (cr) 
starting at xq = 0. Empty circles correspond to a = where 
the numbers label time t = 1, ., 16. Full (and dashed) lines 
represent trajectories for noise amplitude a = 10 -3 (and a — 



io~ 



plotted only at times t = 2" 



as the bifurcation gap 0, Q- For instance, A/j, c ~ a 1 
where Afi = /i c (0) — Mc(c) and 7 = ln<5i?/lnj/, where 
^ = 0.46692... is one of the two Feigenbaum's universal 
constants (the 2nd, a_F = 2.50290..., measures the power- 
law period-doubling spreading of iterate positions), and 
v ~ 2V2a F (l + l/af,)- 1 / 2 ~ 6.619. See @ and refer- 
ences therein. When a is small the trajectories visit se- 
quentially the set of 2 N ^ disjoint bands leading to a cy- 
cle, but the behavior inside each band is irregular. These 
trajectories represent ergodic states as the accessible po- 
sitions have a fractal dimension equal to the dimension of 
phase space. When a = the trajectories correspond to a 
nonergodic state, since as t — > 00 the positions form only 
a Cantor set of fractal dimension dj — 0.538.... Thus 
the removal of the noise a — > leads to an ergodic to 
nonergodic transition in the map. 

As shown in Ref. [f| when fj, = ^ c (cr) (a > 0) 
there is a 'crossover' or 'relaxation' time t x = <J r ~ 1 , 
r = 1 — ln2/lnz/~ 0.6332, between two different time 
evolution regimes. This crossover occurs when the noise 
fluctuations begin suppressing the fine structure of the at- 
tractor as displayed by the superstable orbit with x = 
described above. (See full and dashed lines in Fig. 
For t < t x the fluctuations are smaller than the dis- 
tances between the neighboring subsequence positions of 
the xq — orbit at fi c (0), and the iterate position with 
a > falls within a small band around the a — position 
for that t. The bands for successive times do not over- 
lap. Time evolution follows a subsequence pattern close 



to that in the noiseless case. When t ~ t x the width of the 
noise-generated band reached at time t x = 2 N ^ matches 
the distance between adjacent positions, and this implies 
a cutoff in the progress along the position subsequences. 
At longer times t > t x the orbits no longer trace the pre- 
cise period-doubling structure of the attractor. The iter- 
ates now follow increasingly chaotic trajectories as bands 
merge with time. This is the dynamical image - observed 
along the time evolution for the orbits of a single state 
/i c (cr) - of the static bifurcation gap initially described in 
terms of the variation of the control parameter [i 0. 

The entropy associated to the distribution of the it- 
erate positions within the 2 N bands has the form S c — 
2 N as, where s is the entropy associated to a single band. 
Use of 2 N = t x and t x = a 1 " 1 , r-l ~ -0.3668 @, leads 
to 



t, 



('/Sc 



\(l-r)/r 



(3) 



Since (1 - r)/r ~ 0.5792 then t x -> 00 and S c -» 
as (7 — > 0. See |5| for details on the derivation. We 
have compared this expression with its counterpart in 
structural glass formers, the Adam-Gibbs equation [3|. 



III. 



TWO-STEP RELAXATION 



The time evolution of equilibrium two-time correla- 
tions in supercooled liquids on approach to glass for- 
mation display a two-step process of relaxation. This 
consists of a primary power-law decay in time difference 
t = t% — t\ that leads into a plateau, and at the end of 
this there is a second power law decay that evolves into a 
faster decay that can be fitted by a stretched exponential 
3]. Also, the duration t x of the plateau increases as an 
inverse power law of the difference T — T g > as the tem- 
perature T decreases to the glass transition temperature 
T g . The first and second decays are usually referred to as 
the /3 and the a relaxation processes, respectively 0- An 
observable example of such correlation function, both ex- 
perimentally and numerically, is the Fourier transform of 
the density-density correlation at time difference t. The 
former is known as the intermediate scattering function 
while the latter is known as the van Hove function 

The study of single trajectory properties (one time 
functions) in Ref. y| led to the suggestion that the dy- 
namical behavior in the map at /i c (cr) would show paral- 
lels to the relaxation properties of glass formers. For in- 
stance, the analog of the f3 relaxation would be obtained 
by considering initial conditions xq outside the critical 
attractor since the ensuing orbits display a power-law 
transient as the positions approach asymptotically those 
of the attractor. The intermediate plateau would cor- 
respond to the regime t < t x , described in the previous 
section, when the iterates are confined to nonintersecting 
bands before they reach the bifurcation gap, its duration 
t x grows as an inverse power law of a. The analog of 
the a relaxation was proposed to be the band merging 
crossover process that takes place for t > t x . To explore 
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FIG. 3: Development of time translation invariance (TTI) in 
the correlation c e (ti,t = ti — ti) through the action of noise 
of amplitude a at the onset of chaos. All trajectories start 
xq — and t x is the time to reach the bifurcation gap. 



more closely these similarities we evaluated the two-time 
correlation function 



Ce{t\,t 2 ) = 



(xt 2 x tl ) - (x t2 ) (x tl ) 



, 1 < h < h 



(4) 



for different values of the noise amplitude a. In Eq. Q 
(...) represents an average over an ensemble of trajectories 
all of them starting with initial conditions xq = and 



FIG. 4: Relaxation at the onset of chaos according to the 
correlation function c e (ti,t = t% — ti) defined in Eq. Q for 
an ensemble of 500 trajectories starting at xq = 0. In (a) t\ 
is of the form 2" + 1 and the noise amplitude is a = 10~ 6 
and 10 -7 . Vertical arrows, from left to right, indicate the 
crossover time t x (a). In (b) ti = 2™ — 1 and a — 10 -7 . 



We first address the question of whether the exposure 
of trajectories to noise has the effect of introducing, af- 
ter an initial transient period, a time translation invari- 
ance (TTI) property into the correlation in Eq. Q, i.e. 
c e {ti,ti) — c e (t — t 2 — ii), ti large. The presence of this 
effect would be analogous to thermalization in molecu- 
lar systems after which equilibrium correlations are mea- 



sured or computed (in glass forming systems for T > T g ). 
In Fig. [3] we show how TTI develops and is maintained 
for a sufficiently large time difference interval t = ti — t\. 
The numerical limitation in evaluating accurate values 
for /i c (c) leads to an upper bound for t, but we checked 
that increasing precision in fj, c (cr) leads to a larger in- 
terval for t for which TTI is observed. In view of the 
results shown in Fig. [3|we can conclude that under ex- 
ternal noise of weak amplitude the ensemble of trajec- 
tories 'thermalizes' asymptotically into an 'equilibrium' 
attractor. 

The TTI property in the map still retains certain mem- 
ory of the initial fi, only in a generic way that reflects the 
characteristic symmetries of the period-doubling onset of 
chaos of the logistic map. Notice that in Fig 3 the values 
for ti are of the form t\ = 2" + l. It is found that both the 
initial decay and the value of the intermediate plateau of 
c e (t) are fixed when the values of t\ belong to the se- 
quence t\ = 2" + k, n large with k fixed. However the 
main decay process of c e (t) appears to be independent of 
t%. In Figs. 0]a and^t we show the behavior of c e (t), re- 
spectively, when h = 2 n + 1 and h = 2" - 1, n = 9, 10, ... 
In the first case (Fig. the correlation maintains the 
initial value of unity from t = until the main decay 
process sets in, while in the second case (Fig. 0Jd) there 
is initial decay at short times followed by a plateau that 
ends when the same main decay as in Fig. [3Ji takes place. 
For other values of k the correlation c e (t) shows a behav- 
ior similar to that in Fig. (@Jd) with different values for 
the duration of the initial decay (that we refer to as the 
(3 decay) and for the plateau, but always with the same 
main decay (that we refer to as the a decay). The a 
decay is itself made up of several plateaus the values of 
which alternate when n is varied (as shown in Figs. [3^, 
and QJ)) . In all cases the duration of the main plateau 
coincides approximatively with the crossover time t x at 
which the bifurcation gap is reached. (See the vertical ar- 
rows in Figs. 0Ji and^Jj). Thus, the identification of the 
encounter of the bifurcation gap as the triggering event 
of the a relaxation process Q seems to be confirmed by 
the numerical evaluation of c e (t). 
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IV. AGING 

In glass forming systems when T — T g < (nonequi- 
librium) two-time correlations lose time translation in- 
variance, and the dependence on the two times t\ and 
£2 = t + ti, has the characteristic known as aging 
More specifically, the time scale for the response to an 
external perturbation increases with the waiting time 
t w , the time interval t w = t\ — to from system prepa- 
ration at to = to the moment of the perturbation at t\ . 
As a consequence, the equilibrium fluctuation-dissipation 
relation that relates response and correlation functions 
breaks down Q. In this regime the decay of response 
and correlation functions display a scaling dependence 
on the ratio t/t w 0. 



FIG. 5: Aging according to the correlation c(t w ,t + t w ) given 
by Eq. JHJ for the Feigenbaum attractor (/1 = /i c , a — 0). 
Total observation time is n = 1000. In (a) is shown the 
explicit dependence on the waiting time (from left to right 
t w — 1, 3, 5, 7, 9, 11). In (b) all curves collapse upon 
rescaling t/t m . 



As indicated in Ref. the power-law position subse- 
quences shown in Fig. [21 that constitute the superstable 
orbit of period 2°° within the noiseless attractor at /x c (0) 
imply a built-in aging scaling property for the single-time 
function Xt- These subsequences are relevant for the de- 
scription of trajectories that are at first held at a given 
attractor position for a waiting period of time t w and then 



released to the normal iterative procedure. We chose the 
holding positions to be any of those along the top band 
shown in Fig. [5] with t w = 2k + 1, k = 0,1,.... One 
obtains [j| 



(a) 



xt +tui ^ exp q (-X q t/t w ), 



(5) 



where exp 9 (x) = [1— (q— l)^] 1 / 1-9 and \ q = In a^?/ In 2. 
This property is gradually removed when noise is turned 
on. The presence of a bifurcation gap limits its range of 
validity to total times t w + 1 < t x (a) and so progressively 
disappears as a is increased @. 

When (7 = trajectories are nonergodic and ensemble 
and time averages are not equivalent. For this reason we 
use a time-averaged correlation c(t\,t2) to study aging 
and its related scaling property at the onset of chaos for 
(7 = 0, instead of the ensemble-averaged e e (ti, £2) in Eq. 
(0J. Also for this case c e (t\,t2) is not defined as \u = 
when the initial positions are all xq = or xt m ■ We chose 
for c(£i,t2) the form 



N 



c(t w ,t + t w ) = (1/N) J2 <P {n) (t w )^ n) (t + t w ), (6) 

n—riQ 



where 4>{t) = /^(0) and = 1 — [ix 1 , and with n 

any positive integer and iV > no a large integer. This 
definition of c(t w , t+t w ) is designed to capture the power- 
law patterns of the trajectories at the noiseless onset of 
chaos. Eq. © considers multiples of the two reference 
times t w and t + t w , i.e. times at which trajectories re- 
currently visit a given region of the attractor [l^. In 
Fig. we show c(t w ,t + t w ) for different values of t w , 
and in Fig. 0J> the same data where the rescaled variable 
t/t w = 2™ — 1, t w = 2k + 1, k = 0, 1, has been used. 
We have calculated c(t w ,t + t w ) for different values of 
N and no and found in both cases the same result as in 
Fig. [SJi. The characteristic scaling of aging behavior is 
especially clear. 
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V. SUBDIFFUSION AND ARREST. 

The sharp slow down of dynamics in supercooled liq- 
uids on approach to vitrification is illustrated by the pro- 
gression from normal diffusiveness to subdiffusive behav- 
ior and finally to a halt in the growth of the molecu- 
lar mean square displacement within a small range of 
temperatures or densities ^2 E3 • This deceleration of 
the dynamics is caused by the confinement of any given 
molecule by a 'cage' formed by its neighbors; and it is 
the breakup and rearrangement of the cages which drives 
structural relaxation, letting molecules diffuse through- 
out the system. Evidence indicates that lifetime of the 
cages increases as conditions move toward the glass tran- 
sition, probably because cage rearrangements involve a 
larger number of molecules as the glass transition is ap- 
proached mm. 



FIG. 6: Glassy diffusion in the noise-perturbed logistic map. 
(a) Repeated-cell map (thick dashed line) and trajectory (full 
line), (b) Time evolution of the mean square displacement 
(xf \ — {x t } 2 for an ensemble of 1000 trajectories with initial 
conditions randomly distributed inside [—1, 1]. Curves are 
labeled by the value of the noise amplitude. 



To investigate this aspect of vitrification in the map at 
/x c (cr), we constructed a periodic map via repetition of a 
single (cell) map. This setting has being used to study 
deterministic diffusion in nonlinear maps, in which the 
trajectories migrate into neighboring cells due to chaotic 
motion. For fully chaotic maps diffusion is normal but 
for marginally chaotic maps it is anomalous [T^| . In our 
case we design the map in such a way that diffusion is 
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due only to the random noise term, otherwise motion is 
confined to a single cell. So, we have the periodic map 
xt+i = F(xt), F{l + x) =1 + F(x), I = ... - 1,0, 1, 
where 

F() _ f HcT 2 \+at -l<ar<0, (] 

\1 - n cX 2 \ + a£, 0<x<l. [i) 

Fig. [U< shows the repeated-cell map together with a 
portion of one of its trajectories. As it can be observed, 
the escape from the central cell into any of its neighbors 
occurs when > 1 and this can only happen when 

cr > 0. As a — ► the escape positions are confined to 
values of x increasingly closer to x = 0, and for a = the 
iterate position is trapped within the cell. Likewise for 
any other cell. Fig. HJj shows the mean square displace- 
ment — (xt) as obtained from an ensemble of tra- 
jectories initially distributed within the interval [—1,1] 
for several values of noise amplitude. The progression 
from normal diffusion to subdiffusion and to final arrest 
can be plainly observed as a — ► 0. For small a (< 1CP 2 ) 
\Xf) — (xt) shows a down turn and later an upturn sim- 
ilar to those observed in colloidal glass experiments 01 
and attributed to cage rearrangements. In the map this 
feature reflects cell crossings. 

VI. FINAL REMARKS 

As we have shown, the dynamics of logistic maps at the 
chaos threshold in the presence of noise displays elements 
reminiscent of glassy dynamics as observed in molecular 
glass formers. The limit of vanishing noise amplitude 
(7^0 (the counterpart of the limit T — T g — > in the 
supercooled liquid) leads to loss of ergodicity. This non- 
ergodic state with vanishing Lyapunov coefficient Ai = 
corresponds to the limiting state, a — ► 0, t x — ► oo, of 
a family of small a states with properties reminiscent of 
those in glass formers. Some additional comments may 
be useful to appreciate both similarities and differences 
between the two types of systems. 

Activated dynamics is a standard component in under- 
standing the relaxation mechanisms of glass formation 
0, 0] and so it is pertinent to note that there is a similar 
characteristic in the dynamics that we studied here. It 
has long been known |g, |9( that the addition of external 
noise to a dissipative dynamical system (here a nonlinear 
one-dimensional map) causes its trajectories to escape 
from attractors. For chaotic attractors the mean escape 
time T has an exponential Arrhenius form 

T^T cxp(E /R), 

where E$ and R are the minimum escape 'energy' and 
noise 'temperature', respectively, and To is the inverse of 
the attempt rate. The difference from the usual molecu- 
lar set up is that the escape is from the noiseless attrac- 
tor and not from a minimum in the potential energy. For 
critical attractors, like the onset of chaos, the minimum 



escape energy vanishes and the escape time is expected 
to follow a power law with strong fluctuations 0, llfij . 
We have seen that as a — > an increasing number of 
distinct but small phase space bands are involved in an 
increasingly slow decay of dynamic correlations. On the 
other hand, in the usual picture of activated dynamics 
in glass formers the increase of the relaxation time as 
T — T g — > is thought to be associated to an increase 
in the landscape energy barriers and the assumption of 
some form of cooperative behavior by means of which a 
large number of particles are rearranged through a very 
slow process 

The structure of the phase space regions that are sam- 
pled by trajectories as a function of the noise amplitude 
resemble the manner in which a glass forming system 
samples its energy landscape as a function of tempera- 
ture 0] . The numbers and widths of the attractor bands 
at noise amplitude a can be thought to correspond to 
the numbers and extents of potential energy basins of 
the landscape at a given depth set by the temperature 
T. The 'landscape' for the attractor has certainly a very 
simple structure when compared to that of the multidi- 
mensional energy landscape in a molecular system. At 
(7 = there is an infinite set of minima that consists of 
all the points M — > oo of the attractor, and as a increases 
these minima merge into finite sets of M = 2 N ^ bands 
with N decreasing as a grows. The regular merging by 
two in the numbers M and the features in the dynam- 
ics that they imprint are properties specific of logistic 
type maps. Other kinds of critical multifractal attractors, 
such as those for the critical circle map 6] , would exhibit 
other properties characteristic of the route to chaos in- 
volved. This availability of detail would not be generally 
present in the measurements or numerical computation 
of the dynamics of a glass forming molecular system. 

It is of interest to note that at fj, c (a) the trajectories 
and its resultant sensitivity to initial conditions are ex- 
pressed for t < t x via the g-exponentials of the Tsallis 
q-statistics [5j- For a = this analytical forms are exact 
fiol Il5j and an identity linking accordingly generalized 
Lyapunov coefficients and rates of entropy production, 
holds rigorously [ill Il5| . There is nonuniform conver- 
gence related to the limits a — > and t — > oo. If a — > 
is taken before t — > oo orbits originating within the at- 
tractor remain there and exhibit fully-developed aging 
properties, whereas if t — > oo is taken before a — > a 
chaotic orbit with exponential sensitivity to initial con- 
ditions would be observed. 
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